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ABSTRACT 

The model of non-ballistic jet motion proposed in 2008 provides a simple explanation 
to the inward jet motion and bent jet. Recently, evidences of such a non-radial motion 
increase rapidly, and more complicated morphologies appear. On the other hand, the 
ballistic plus precession model likely holds in majority samples of jet motion. This 
paper discusses the relationship between the ballistic and non-ballistic model of jet 
motion, which suggests that the interaction of ejectors with ambient matter can pro- 
duce knots at different stages of evolution and hence different separations to the core. 
And as a jet precesses, knots produced between the core and the deceleration radius 
result in spiral pattern expected by the model of ballistic plus precession; and knots 
generated at the deceleration radius display non-radial motion such as bent jet or os- 
cillation of ridge-line. This paper develops the first non-ballistic model in four aspects. 
Firstly, it provides a numerical simulation to the production of multi-knot for a pre- 
cessing jet. Secondly, it fits the precession behavior of multi-knot and interprets the 
oscillation of ridge lines like S5 1803-1-784. Thirdly, it gives an unified interpretation to 
the bent jet applicable to both multi-knot and single knot. And fourthly, the problem 
of very large numbers of observed outward motions as opposed to the inward ones is 
addressed in a new scope. 
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1 INTRODUCTION 



To explain phenomena of non-radial jet motion, which are 
difficult to understand under the context of ballistic motion 
of out flow ( . Ree^ 1966 ). the non-ballistic model was pro- 
posed (|Gonel |2008| ) . in which a continuous jet produces a 
discrete hot spot by the interaction of jet-matter at the de- 
celeration radius. And the precession of such a hot spot in 
the plane of the sky interprets the inward and bent motion 
of AGNs (|Kellermann et al.ll2004l : lAgudo et al.]|2007l ). 

Recently, more and more outflows showing non-radial 
motion were observed. For example, the MOJAVE program, 
which is intended to investigate the parsec-scale jet kinemat- 
ics of a complete flux-density-limited sample of 135 active 
galactic nuclei (AGN), has revealed many unusual proper- 
ties of AGN jets. Interestingly, about one third of the outflow 
components have velocit y vectors tha t do not point ba ck to 
the core feature (jHoman et al. 2009: . Lister et al]|2009l ). 

Complicated non-radial phenomena were also ob- 
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served in a set of interest ing BL Lac Objects e.g., 
BL Lac (IStirling et all |2003|) 0716 -H714 (|Britzen et all 
l2009l) . 1803 +784 (iBritzen et al] 2010al ) (1803 hereafter) and 
0735-1-178 (|Britzen et all l2010bl ). which display oscillation 
behavior of ridge lines. While it is usually expected that the 
outflow components should escape from the central core at 
apparent superluminal speeds of about 5 — 30c, the new 
observations actually show that in about 20 years or longer, 
the knots appear stationary with respect to the core while 
exhibiting signiflcant variation in their position angles. As a 
result, the ridge lines of 1803 change substantially. Such an 
oscillation of jet ridge line associating with the production 
of multi-knot in the jet-matter int eraction has not been ad- 
dressed in the non-ballistic model (|Gonell2008l ). which needs 
to be developed. 

Interestingly, the non-radial motion displayed in AGNs 
is also observed in the famous X-ray binary system of SS 433, 
where a pair of reversely moving, mildly relativistic jets are 
detected at X-ray, optical and radio wavelengths l|Margonl 
[l984). In t he 2001 observation campaign on the arcsec-scale 
X-ray iets (|Migliari et al.ll2005l ). a knot (knot A) brightened 
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from May 8 to May 10, and then seemed to be moving along 
a precession trace to a new place (knot B) on May 12. Again, 
the behavior of these two knots is not easy to understand 
in the framework of the standard ballistic model. If both 
knot A and B are produced by a single outflow traveling 
from A to B, then the position angle should not change so 
dramatically as observed, and the travel time also should be 
much longer than the reported 2 d. On the other hand, if 
they are produced by two discrete outflows, then B should be 
ejected earUer than A, and should correspondingly appear 
earlier. This is contrary to the observations. Moreover, the 
ballistic scenario predicts that both knot A and B should 
move outward as a whole, which is also not observed. 

In some cases the jet motion like knots A and B of 
SS433, have been used as evidences contradicting to the 
scenario of jet precession. In fact, the problem is originated 
from the ballistic assumption rather than the jet precession. 
Therefore, the relationship between non-radial and radial jet 
motion needs to be discussed. 

Non-radial jet motions are difficul t to underst and in the 
context of ballistic motion of outflows |Rees|[l9 66l). but they 
can be reasonably explained by the non-ballistic model as 
proposed by Gong (2008), who argued that a continuous jet 
can produce discrete hot spots at the deceleration radius due 
to its interaction with the ambient matter. The precession of 
the jet can then lead to apparent motion of the hot spots on 
the sky plane, which explains the inward and bent motions 
of some AGN knots. 

In this paper, we go further to show that the interac- 
tion between the jet and ambient matter can give birth to a 
number of knots along the outflow. The knots can be gen- 
erated either at the deceleration radius or much closer to 
the core. Knots produced near the deceleration radius can 
naturally ex plain the variation of ridge line s as observed 
in BL Lac (Istirhng et al.1 120031 ) and 1803 l|Britzen et al.1 
l2010al ). On the other hand, knots produced closer to the 
core can explain the spiral pattern when the precession of 
the jet is considered. The structure of our paper is as fol- 
lows. In Section 2, the difference and relation between the 
shock-in-jet model (Marscher & Gear 1985) and the non- 
ballistic model is addressed. The mechanism of producing 
multi-knots is described in detail. In Section 3, the model 
is engaged t o explain the oscilla tion of ridge lines of 1803 
and BL Lac (|Stirling et al.ll2003l ). as well as the non-radial 
motion of SS 433 knots. In Section 4, we discuss the com- 
mon mechanism underlying these sources. A simple formula 
is derived to describe the curvature of non-radial motions, 
which provides a general insight int o the non-radial behavior 
displayed in various AGN sources (|Kellermann et al] l2004l : 
iHoman et aLll2009l : iLister et al.ll2009l '). Finallv. Section 5 fo- 
cus on the following interesting problem: why the number 
of outward-moving knots are much larger than that of the 
inward-moving ones. 



2 THE NON-BALLISTIC MODEL IN 
MULTI-KNOT CASE 

2.1 Shock-in-jet model and non-ballistic model 

The flux density of a knot (knot i) can be expressed as 

s,{u,R„e) = di^D{R,,ef~^j[{u)vi (1) 



where 6 is the viewing angle between the jet axis and our 
line of sight (LOS), which causes variation in the flux den- 
sity due to the Doppler boosting effect. The Doppler factor 
is given by D{Ri,9) = 'y"^ {Ri)[l - (]b{Ri) cos 9]-'^ , with the 
Lorentz factor 7(,(i?,) = [1 - l3'i{Ri)]'^^'^ . Vi in Eq. (1) is 
the volume in the co-moving frame of knot i, which is con- 
nected with the volume in the observer's frame by Vi = DV^ . 
ji is the emissivity, which can be derived by considering the 
detailed mechanism. There are mainly three kinds of mod- 
els to explain the outflow radiation, the shock-in-jet model 
(Marscher & Gear 1985), the internal shock model (Spada 
et al. 2001), and the non-ballistic model (Gong 2008). We 
discuss them one by one below^ 

The shock-in-jet model (iMarscher fc Gearl Il985l ) ex- 
plains the light-curve (L-C) by assuming that each peak in 
the L-C is associated with the evolution of a single shock 
front. In this model, the shock wave propagates along the 
jet, experiencing three different regimes sequentially, i.e. 
the inverse Compton loss regime, the synchrontron radia- 
tion regime, and finally the adiabatic regime. At each stage, 
the emissivity, , can be calculated when the parameters of 
radius, Ri, bulk velocity, /3(,(-Ri), matter density, magnetic 
field, etc, are known. 

On the other hand, the internal shock model of gamma- 
ray bursts (GRBs) has been ge neralized to accou nt for the 
outfiow evolution. In this model ( Spada et al ]|200ll ). the cen- 
tral engine producing the outflows is assumed to work in- 
termittently, so that later but faster shells can catch up and 
collide with slower earlier ones, producing strong shocks. 
The shocks can convert a part of the bulk kinetic energy 
into internal energy of electrons which then can be dissi- 
pated via various radiation mechanism to manifest as bright 
knots. 

Note that the knots produced in both the shock-in-jet 
model and the internal shock model should be very close to 
the core, at least much closer as compared with the deceler- 
ation radius. However, the deceleration radius, used by the 
non-ballistic model, is a natural consequence of the dynami- 
cal interaction between the outflow and the circumambience 
matter. It is very similar to the external shock model widely 
used in investigating the afterglows of GRBs. Assuming that 
a beamed outflow with a kinetic energy E is ejected from 
the core. It expands into the ambient medium with a den- 
sity of no. The supersonic motion of the ejecta should drive a 
blast wave propagating into the interstellar medium (ISM). 
At the same time, the ejecta itself should be decelerated due 
to more and more swept-up material. 

The evolution of the outflow can be calculated by 
considering the co nservation of energy (IWang et al.l l2003l : 
iHuang et ahlbOOOl '). which can be simplified as, 

Eshocki (tl) + Eshelli (tl) = E{t\) . (2) 

where ti is the starting time of the activity, E{t\) and 
Eshockiiti) are the total energy of the outfiow and the ki- 
netic energy of the shock, respectively, and Esheiiiiti) is 
the internal energy of the shock. The deceleration radius is 
defined as the radius where the mass of the swept-up mat- 
ter equals to 1/71, of the initial mass of the outflow. The 
external shock, and consequently the emissivity, should be 
the strongest at the deceleration radius. So it also deter- 
mines the positions of the observed knots. Therefore, in the 
non-ballistic scenario, the bulk velocity, l3b{Ri), Ri, and ji 
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of Eq. ((TJ, can be obtained from the dynamics described 
by Eq. ((2}. Note that the deceleration radius is mainly de- 
termined by the energetics of the outflow and the density 
of the ambient medium. It means that the distance between 
the knots and the central core should also be relatively fixed. 

The shock-in-jet model can produce a sequence of knots 
by adopting repeated onsets of several shocks. Consequently, 
when the jet precesses, a spiral pattern consisted of the knots 
may be formed. It is very similar to the ballistic plus preces- 
sion scenario previously used to explain the morphology of 
SS 433. In other words, the shock-in-jet model itself is one 
kind of ballistic plus precession model, owing to the bulk 
speed of the knots. 

In contrast, in the non-ballistic model, a knot should 
stay at approximately cons tant separat ion (at the decelera- 
tion radius) from the core (|Gongj|2008l ). When the jet pre- 
cesses, a sequence of knots can also be produced, but all the 
knots should be static. They should not move significantly 
outward. Of course, since the matter density at different 
directions may differ, the separations of the knots from the 
core is not necessarily a constant at different directions. This 
provides a possibility of varying the input parameters when 
we calculate the fiux of a knot by using Eq. Q. 

In short, the non-ballistic model differs from the shock- 
in-jet model in at least two aspects: (i) It can produce knots 
that reside at deceleration radius, (ii) It can naturally gen- 
erate a ring-like trace, instead of a helical one as predicted 
by the shock-in-jet model. 

2.2 The production of multiple knots during jet 
precession 

Eq. ^ is applicable for the continuous jet of AGNs and 
X-ray binaries thanks to the precession of the jet. The ac- 
tive timescale of the central engine is usually negligible as 
compared with the precession period. For example, for a 
precession period of 10 years, the typical burst timescale of 
St — Iks is sufficiently short, during which the change of 
the jet precession phase is only ~ 3 x 10~^rad. In such 
a short time scale, the energy ejection can be regarded as 
instantaneous and the interaction with ambient matter can 
be treated by using Eq. (Q. It will appear in the sky as a 
bright knot. 

At another moment of t2 = ti + 5t, another ejecta comes 
out at a slightly different direction with respect to the LOS, 
02 = 9i + 50. Such a process can occur repeatedly, giving 
birth to multiple knots. 

If the jet is continuous and with exactly constant power, 
E{ti) =const, and the surrounding ISM is homogeneously 
distributed, then the process described by Eq. ([2} may per- 
form numerous times (e.g., 10yr/lks=3 x 10^), and knots 
should be produced farther and farther away from the core. 
However, many AGNs show variability of t ime scale of 
- lOks (|Gliozzi et al.l 12004 [Rani et alj|2010l ). The begin- 
ning of such a pulse, with a period of ~ lOks, means restart- 
ing the process of Eq. ((2]) again. 

Consequently, during a pulse (not necessarily a con- 
stant), a numbers of bursts, say, 10, are generated, and 
only a few of them can reproduce L-Cs as shown in Fig. 1. 
Other bursts may produce knots at places deviating sig- 
nificantly from the deceleration radius or with flux density 
much weaker than that in Fig. 1. Thus, statistically the pat- 



tern formed at each pulse is stable, in which the knot-core 
distance of a knot is approximately constant or constrained 
in a certain range. A schematic demonstration of such knot 
production is shown in Fig. 2, in which X1-X3 correspond to 
knots produced in one pulse, and Y1-Y3 by another pulse. 
The upper panel differs from the bottom one in precession 
phase, the former has small discrepancy and the latter one 
large. 

Therefore, at a direction slightly differs from the previ- 
ous one, once the energy ejection of the pulse and the matter 
density etc satisfy the parameters of Table 1, knots with L- 
Cs as shown in Fig. 1 can be reproduced again. The small 
phase discrepancy among the knots produced in one pulse of 
the central engine results in the pattern like, knots X1-X3, 
and large phase discrepancy corresponds to the pattern of 
Y1-Y3 in Fig. 2. 

On longer time scale, e.g., with 10^ pulses, the change 
of jet precession phase is of ~ lOdeg, during which the mat- 
ter density can be variable through e.g., evacuated bubbles 
around the sources due to previous activities of the jets. By 
Eq. ((2)1 such a fluctuation in matter density makes the L-C 
peak at different distance to the core as shown in Fig. 1, 
which varies the knot-core separation of knots at different 
directions. Thus, the ring through XI and Yl, as shown in 
Fig. 2, can be tilted, so are rings through X2,Y2 and X3,Y3. 

Therefore, the formation of multi-knot morphology can 
be realized by extending the former non-ballistic model, 
where several knots are produced along the jet axis. And by 
the jet precession, these knots cause ring like traces at dif- 
ferent separation to the core. Observing at different epochs 
result in the oscillation of ridge lines. 

The ring like pattern corresponds to the knot produced 
at the deceleration radius. Comparatively, spiral patterns 
correspond to the knots produced between the deceleration 
radius and the core, which can be interpreted by the shock- 
in-jet model. 

Thus simply replacing the knot-core distance, Rd of 
Gong(2008), describing the precession of one knot under the 
non-ballistic model by the multi-knot, (where i = 1, 2, 3... 
corresponds to different knots). The projection of Ri to AS 
and Aa axes gives, 

Ra = R^ [sin A sin r)^] sin ^ + R^ [cos A sin I — sin A cos / cos r\ \ cos ^ , 
-Ri = K [sin A sin r\ \ cos — Rl [cos A sin 7 — sin A cos / cos r] \ sin($) 

where A is the opening angle of the precession cone, 7 is the 
inclination angle between the jet rotation axis and LOS. The 
precession phase is jf — Clt + rjo {fl is the precession velocity 
of the jet, rjo is the initial phase of each knot). Projecting 
R^ into the coordinate system x ~ y ~ z, where the a;-axis is 
towards the observer. Then rotating around the a;-axis for 
angle ^, so that the new y-axis (AS) will point north, and 
the new z-axis (Aa) will point east. The position angle of a 
knot can be simply obtained by 



= tan~\Aa/A5) . (4) 
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3 APPLICATION TO AGN AND X-RAY 
BINARY 

3.1 Oscillation of ridge line 

The blazar S5 18 03+784 is a flat-spectrum radio source at 
high declination (|Witze]||l987l ). Geodetic and astronomical 
VLBI data gathered at 8.4 and 5 GHz between 1979-1987 
showed that the co mponent located at 1.4 mas from th e core 
appears stationary (jSchalinski et al.ll2003l : [Witzel 19881 ). The 
stationary component was found to ha ve non-constant cor e 
separation or oscillatory type behavior (|Britzen et al.|[2005h . 
And rece ntly the oscillation of ridge lines of S5 1803-1-784 is 
reported (|Britzen et al.ll2010al '). 

In the following we describe the fitting of the kinemat- 
ics in the pc-scale jet of 1803 by the extended non-ballistic 
superluminal model of Section 2. The evolution of core sepa- 
ration, flux and position angle of one compon ent, CI, in the 
time in terval between 1984-1996 are observed l|Britzen et al.l 
l2010al ). as denoted by the flUed squares shown in Fig. 3. Un- 
der the non-ballistic scenario, the behavior of CI is simply 
a conical precession of one knot, projecting to the plane of 
sky. 

With the fitting parameters of Table 1, the peaks of L- 
Cs in Fig. 1 can accou nt for the observatio nal knots, e.g., 
Cl,Ca and C4 of 1803 l|Britzen et al.ll2010al ). By Fig. 1, at 
a time scale of half a year, the peaks can decline for 20%- 
10%. More rapid declination is expected in case of larger 
spectral i ndex of electron , po, or a radiative dynamics of 
the shockl Sari et al.l Il998l ). And considering the effect of 
Doppler boosting, the L-C of a knot can vary more dra- 
matically. 

The fiux fitting of Fig. 3 corresponds to a bulk speed 
of Pb{Ri) ~ 0.3 — 0.6. Therefore, the core separation, posi- 
tion angle, and flux density of CI, can be fitted by Eq. Q, 
Eq. Q and Eq. ((T)) respectively. The variation of sep ara- 
tion, position angle and flux of CI l|Britzen et aLllioiOal ) are 
fltted by two groups of free parameters. Group A contains 
9 global parameters, such as precession speed of jet axis, 
the opening angle of precession cone and the bulk speed as 
shown in the flrst row of Table 2. Group B includes 13 os- 
cillation parameters, denoting the deviation of the distance 
of component CI to the core at 13 different epochs from the 
averaged value as shown in Table 2. These two groups of 
parameters totally 22, fit three figures, the evolution of the 
core separation, flux and position angle at 13 epochs, totally 
13 X 3 = 39 observational points. The three fltting results 
are reasonable well, as shown in Fig. 3. Notice that in the 
fitting of position angle and flux of Fig. 3 (26 observational 
points) only the 9 global parameters are needed, as shown 
in the first row of Table 2. 

Other components, C2,C3...can be treated similarly, ex- 
cept discrepancies on the knot-core separations and phases. 
The projection of these knots at di fferent time explain s 
the oscillation of ridge lines observed (|Britzen et al.|[2010al ). 
Variation of knot-core distance up to 50 percent of the aver- 
age one, Ro, is needed in the fitting of Fig. 3, which indicates 
the fluctuation of ISM distribution and hence large discrep- 
ancy in knot-core distances at different directions. 

The time taken by a knot to precess for tangent distance 
of the size of a knot, Sr, is Stpr = Sr/r/Cl, where Srjr = 
5r/{R^ sin A). The cooling time of a knot, Stco, corresponding 
to, i.e., the time taken for the radio peak to decline for a 



certain percentage, can be inferred from the L-C, as shown 
in Fig. 3. 

3.2 Precessing jet nozzle 

Very Long Baseline Array (VLBA) in an observing mode 
sensitive to linear polarization at wavelength 7 mm with a 
resolution of order of 0.2mas has been perfor med on BL Lac 
at 17 regular epochs from 1998.23 to 2001.28 l|Stirling et all 

liooay 

The observations suggest relatively straight trajectories 
near the core, increasing in curvature at large separations to 
the core(greater than 1-2 mas). The observed continuation 
of the trajectory of component SIO did not fit the prediction 
of helical model jPenn et al.l|2000| ). 

As shown in Fig. 4, ridge lines formed by four knots 
are obs erved in Decem ber 1999, April 2000 and April 
2001 (S tirling et ai]|2003l ). which are labeled by points with 
error bars, A1-A4, B1-B4 and C1-C4 respectively. 

The largest discrepancy appears in the ridge line in De- 
cember 1999, when the helical model predicts a ridge line 
approximately along the vertical line from (0.0, 0.0) to (0.0, 
-0.6), the observed ridge line corresponds to points with er- 
ror bars, A1-A4 in Fig. 4. 

By the non-ballistic model, with fitting parameters of 
Table 3, the observed knots, the ridge lines A1-A4, B1-B4, 
and C1-C4, can be well fitted by the filled circles, squares 
and triangles respectively, as shown in Fig. 4. The relatively 
small discrepancies in a few points in Fig. 4 can be im- 
proved by assuming small variation in the core-knot sepa- 
ration. This means that even the straight ridge lines can be 
explained by knots predicted by the non-ballistic scenario. 
Consequently, the true ballistic jet motion should occur in 
the region even closer to the core than these straight ridge 
lines of Fig. 4. 

It was suggested that initial straight component trajec- 
tories and the subsequent bending jet is due t o a transition 
from a ballistic fashion to non-ballistic flow (|Stirling et al.l 
I2003D . 

Whereas, from the point of view of non-ballistic model, 
the non-radial knot can be interpreted simply by extending 
the fitting of Fig. 4 to knots with larger separation to the 
core and with larger discrepancy in the initial precession of 
each knot, which is further discussed in Section 4. 

3.3 Two puzzling knots of SS 433 

Interestingly, the non-radial motion displayed in AGNs is 
also observed in the famous X-ray binary system of SS 433. 
, in which anti-parallel jets traveling near-relativistically 
are detected in observations across X-ray, optical and ra- 
dio wavelengths (Margon, [1984). The strong and broad 
emission lines have been identified as redshifted/blueshifted 
emission from coUimated jets with intrinsic velocity of 
~ 0.26c jAbellfc MargonI 1 19791 '). The periodic change of 
Doppler shift of emission lines with time is widely accepted 
to be a precessing ballistic jet with a period of 164d. 

Although the jet motion of SS 433 is less than the speed 
of light, the behavior of which is strikingly similar to the 
superluminal sources. In the 2001 observa tion of arcsec-scale 
X-ray jets of SS 433 (|Migliari et al.ll2005l 'l. As shown in the 
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top panel of Fig. 5, a knot became brighter from May 8 to 
May 10 (knot A) and subsequently on May 12 it appeared 
to be moving along a precession trace to a new place (knot 
B). 

The observation of the two knots seems simple, it is not 
easy to understand in the context of the standard ballistic 
plus precession model either. If knots A and B are results 
of one ball traveling from A to B, then its position angle 
should not change so dramatically, and the discrepancy in 
the travel time of these two knots to the observer should be 
much longer than the reported 2 d. On the other hand, if 
knots A and B are caused by two ejections, then B should be 
ejected earlier than A. As a result, B should appear earlier 
than A. This is contrary to the observation. Moreover, the 
ballistic scenario predicts that both knots A and B should 
move outwardly as a whole. Whereas, one cannot find such 
a motion either. 

To explain the motion of these two knots under the stan- 
dard ballistic model, the underlying faster outflow scenario 
is proposed, in which the two knots ejected from the binary 
core first, and later a fast shock wave propagates from the 
cor e through the jet an d hits first the knot A and then knot 
B (|Mid lari et al.ll2005l ). This requires the central engine to 
switch on different modes of power at different time. 

In some cases the jet motion like knots A and B of 
SS433, have been used as evidences contradicting to the 
scenario of jet precession. In fact, the problem is originated 
from the ballistic assumption rather than the jet precession. 
Therefore, the relationship between non-radial and radial jet 
motion needs to be discussed. 

The model of two knots moving at certain posi tion and 
then brightened up by two shocks sequently ( Migliari et al.l 
[iooi), is shown in the middle panel of Fig. 5. Knot B is at 
a distance of 1.7 x lO^^cm from the core, which was ejected 
from the binary system about 40d earlier than knot A, and 
is moving towards us with an angle to LOS of ~ 80deg. The 
shock wave interact first with knot A and then with knot B. 
When the shock interacts with knot A, the distance between 
A and B is about 20 light days. Therefore, the shock wave 
has to travel the projected distance between the two knots 
in about 22d. In order to observe the brightening of the two 
components within two days, the shock has to travel with a 
velocity oi v ^ 0.5c. 

The non-ballistic model provides a simpler scenario, 
in which knot A with a X -ray peak count of 75.2 at ~ 
10 cm (|Migliari et al.ll2005l ) can be reproduced at the de- 
celeration radius (measured from the central), where the 
swept-up medium b y the jet has an energy c omparable to 
that of the outflow jBlandford fc MckedflQT^ '). 

rd = {- — 2 7)^ (5) 

where 70 is the initial Lorentz factor of the outflow. So the 
separation of knot A to the core of ~ lO^^cm corresponds 



to ^5^70,5^' 



0.54, where E51 is E in units of 10^^ erg. 



70,5 is 70 in units of 5, no is in units of 1 cm~ . 

By the non-ballistic model, knots A and B observed in 
SS 433 can be fltted by the precession trace as shown in 
the bottom panel of Fig. 5. The fitting parameters are given 
in Table 2. Therefore, knots A and B can be lying on the 
precession trace given by Eq. Q. The deviation in precession 
phase of A and B corresponding to a time discrepancy of 38 



days. It takes 36 days for a signal to propagate from A to 
B. Hence knot B is brightened up two days later than knot 
A, which satisfies the observation (jMidiari et al.ll2005h m a 
much simpler way. 

The fitting of the two dimension trajectory of knot A 
and B by Eq. Q actually needs 4 equations (corresponding 
to the coordinate of A and B), and the 2d deviation in prop- 
agation time imposes another constraint on the discrepancy 
of the knot-observer distance. Therefore, the bottom panel 
of Fig. 5 can be recognized as obtaining the five variables, 
as shown in Table 4 by five equations (constraints). 

Note that in principle, no matter how many input pa- 
rameters are used, the trajectory of knot A and B cannot 
be fitted by the ballistic plus precession scenario directly. 



4 THE NATURE OF BENT JET 

The origins of jet curvature in the one-knot case has been 
discussed in Gong (2008). Here it is extended to the multi- 
knot case, and the association of the two types of bent jets 
is addressed. 

Within the non-ballistic model, there are two origins 
of jet curvature. The first one is caused by the knots mov- 
ing at different rings. In this case, the phase discrepancies 
between two neighboring knots at different separations are 
given, Srj^ — tlSt^ + in which Sf^ can be written. 



5f 



t — t — oR /v + [x 



x']/c, 



(6) 



where x'' is the knot-observer distance of knot i, 5R^ = 
R^^^ — R^ is discrepancy of the core-knot distance between 
knots i and i + 1, and is the bulk speed of knot i -I- 1. 
A dramatic reduction of bulk velocity between e.g., knot 
2 and 3, « , results in SR^/v^ > SR^/v^, which 
leads to a large phase discrepancy, Sr/^ between knot 2 
and 3. In such case, the jet is bent sharply at knot 2 by 
Eq. ((6]) and Eq. JS]). This explains the wiggling structure 
of 0548-1-165 (|Mantovani et al.iriQQSl '). and the two rectangu- 
lar pa,t tern in BL Lac object PKS 0735-1-178 (|Britzen et all 
l2010bl ). 

Knots produced in the case of small phase discrepancy 
and large phase discrepancy can be demonstrated by the 
upper and bottom pattern of Fig. 2 respectiv ely, which ex- 
plains the variation of ridge lines of BL Lac (|Stirling et al.l 
|200a) and the oscillation of ridge lines of 0735-1-178 and 1803 
respectively. 

The second origin of jet curvature corresponds to the 
trajectory produced by a singular knot, with negligible vari- 
ation in radial separation, which is simply part of an ellipse 
in th e plane of the sky, which has been discussed (|Gonel 
120081 ). 

From the stand of Eq. ((6|, the curvature in such one- 
knot case can be obtained simply by removing the first term 
at the right hand side of Eq. ((SJ, which well explains the 
curvature of knots A and B in SS 433. 

In fact, such a singular knot curvature determined by 
Eq. © and Eq. ^ can explain the non-radial behavior re- 
vealed in other AGN sources, where outfiow are aligned with 
the local jet direction, suggesting the jet flow occurs along 
pree xisting bent channels, like 0738-1-313 and quasar NRAO 
150 jKellermann et al.|[20oi : iLister et al.ll2009l : lAgudo et af] 
1200711 
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By the time derivatives of Eq. (|3]), and R} can 
be obtained, which naturally provides the acceleration per- 
pendic ular to the velocit y vector displayed in AGN compo- 
nents (|Lister et al.ll2009l '). 



5 WHY 98.5% OUTWARD MOTION 

Although the problem that the jet features showing an over- 
whelming tendency of outward motion away from the core 
has been addressed (|Gonell2008l ). it can be further discussed 
in three aspects. 

The first is that most jets of AGNs and X-binaries are 
observed at the stage before ejectors reaching the decelera- 
tion radius, which is dominated by the ballistic scenario. 

The second is that some sources have a small opening 
angle of precession cone, the geomet ry of which displays the 
oscillation of ridge line like BL Lac (jStirling et al.ll20C)3r ). In 
such case, outward and inward motions are performed in 
limited region, and motion perpendicular to the ridge line 
is obvious, which is usually cataloged to the helical motion 
instead of inward motion. 

The third one is for those sources that the jet axis really 
through LOS during the precession of jet. In such case, fol- 
lowing mechanism can make asymmetry in the inward and 
outward motion of such a precession jet. 

By Eq. ((5]), a mildly curved jet appears when SR^ /v'^^ 
is not too large at different locations of the jet. As shown in 
Fig. 6, a knot is produced at the deceleration radius at a pre- 
cessing and curved jet. Since the knot actually corresponds 
to a region, e.g., the knot CI of 1803 is a bout 5r ~ 0.5 
mas(3pc) in diameter (|Britzen et al.ll2010al ). Owing to the 
energy dissipation at the deceleration radius, the bulk speed 
of ejector reduces to typically half of its original speed, which 
can be simplified as composed of several sub-knots, each cor- 
responds to a sub-jet as a and b in the left panel of Fig. 6. 
According to the standard GRB model, a dramatic reduc- 
tion of bulk speed occurs at the deceleration where the knot 
produced (e.g., from 7 = 10 to 7 = 2), so that the upstream 
sub-knots posses larger bulk speed than those of the down- 
stream ones, 7i > "/i+i, by a factor of a few. Hence through 
9i = 1/7j, the half opening angle of emission beam of these 
two sub-knots, i and i + 1, satisfy, 9i < 9i+i. Thus the ef- 
fective emission beam of the knot is extended, with higher 
brightness at upstream region and lower at the downstream 
region, as shown in the middle of Fig. 6. 

With a represents the misalignment angle between two 
sub-jets a and b corresponding to the two knots having the 
minimum and maximum beam opening angles, dmin and 
e^a^. E.g., for knot CI of 1803, a ^ 5r/Ro = 0.5/0.8 = 
0.6rad. As shown in the left panel of Fig. 6, a knot beam 
extended significantly when a ^ 9m,in{e.g.,9min ~ 0.1, with 
7 Ri 10); and a negligible extension occurs when a « 9min 
{e.g.flmin « 0.5, with 7 « 2). 

The precession of the curved jet of the left panel of 
Fig. 6, is equivalent to a fix jet and a change of the observer's 
LOS, as shown in the middle of Fig. 6. At epoch 1, the ob- 
server sees a brief inward motion, and epochs 2-4 correspond 
to the long outward motion due to the extended emission 
beam, as shown in right panel of Fig. 6. Thus the very large 
numbers of observed outward motions as opposed the in- 
ward ones results when most of the knot emission beams 



are deformed considerably. Consequently this causes serious 
asymmetry at two sides of the line connecting the observer, 
the knot and the core (case 2 in the middle and right hand 
side of Fig. 6) . 

Notice that according to this scenario, the size of an 
inwardly moving knot is not necessarily large than that of 
an outward one. Whereas, the explan ation of inward motion 
by jet bending back and across LOS fMars cher et al.lll99ll ) 
predicts larger size of inward knot than that of outward ones, 
considering the lateral expansion of a knot at very longer 
distance from the core. 

In this paper, multi-knot corresponds to a series of knots 
produced by jet-matter interaction. Singular knot is an in- 
dividual knot produced by the jet-matter interaction. And 
sub-knot represents a singular knot (or one of the multiple 
knots) at different stage of evolution near the deceleration 
radius. 



6 DISCUSSION AND CONCLUSION 

This paper shows that generally the interaction of jet and 
ambient matter can reproduce a number of knots along 
the jet, both at the region near the deceleration radius 
and at the region between the deceleration radius and the 
core. Knots produced at the region near the deceleration 
radius are simulated by the dynamics of jet-matter inter- 
action. And the motion of such multi-knot is used to ex- 
plain the changes of ridge line such as BL Lac (|Stirling et al.l 
I2OO3I ) and 1803 l|Britzen et aLlbOlOal ). On the other hand, 
knots produced between the deceleration radius and the 
core ca n be understood in the context of the Shock-in-jet 
model (|Marscher fc ^GeaBIHii), in which a spiral pattern 
is expected when the jet precesses. In other words, this pa- 
per not only extends the former non-ballistic model from 
one-knot to multi-knot, but also discusses the relationship 
between the non-ballistic model and the Shock-in-jet model. 

Interestingly, jets observed by the early long base line 
technique e.g., in 1970s-1980s, appear more straight than 
those observed in the past ten years with much improved 
sensitivity. From the view of non-ballistic model, this is ex- 
pected. Because even a curved jet can behave as a straight 
one, due to the low sensitivity observation favors to observe 
bright knots having small misalignment between the jet axis 
and LOS, which is Doppler boosted strongly. In other words, 
it is this selection effect that makes the early jets appear 
straight. By the improved sensitivity, weaker and weaker 
part of jet is measurable, so that more and more curved jet 
appear. Thus, the ratio of non-radial over radial motion of 
jet should increase with time. 
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Table 1. Tlic fitting parameters of three knots by the standard 
afterglow of GRB model. 



knot 


To 








Eo 


no 


Po 


1 


5.0 


0.1 


0.01 


0.01 


6.0 X lO^f 


0.2 


2.1 


2 


5.0 


0.1 


0.01 


0.02 


2.0 X 10^^ 


0.02 


2.3 


3 


5.0 


0.1 


0.01 


0.03 


8.0 X lO^s 


0.002 


2.8 



The initial isotropic equivalent kinetic energy of the outflow, Eq 
(erg), opening angle of jet, 9j (rad), initial Lorentz factor, Fq, 
density of ambient matter, n, spectral index of electron, poi the 
fraction of total shock energy acquired by the shocked electron 
and magnetic field, Ee and eb respectively. 



Table 2. The fitting parameters of the knot-core separation, po- 
sition angle and flux density of component CI of 1803 as shown 
in Fig. 6. 
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A 




Ro 


/36 


n 


50 


1.43 


36.3 


1.27 


0.28 


4.24 


0.82 


0.34 


2.00 


0.10 


5Ri 


SR2 


SRs 


5R4 


SRb 


SR6 


SRr 


SRs 


SR9 


-0.37 


-0.45 


0.12 


0.17 


-0.018 


-0.16 


-0.31 


-0.20 


-0.05 


SRw 


SRii 


SRi2 


SRia 












0.00 


-0.046 


0.060 


-0.37 













The average knot-core separation, Ho, and the deviation of the 
knot-core distance from the average value at different time, SRi 
(i = 1 — 13), are in unite of mas; Q is in deg/yr, A, I and rjo 
are in rad. The two parameters, n and So denote the index 3 — a 
and the flux density at the righthand side of Eq. Q respectively. 



Table 3. The fitting parameters of BL Lac (Stirling et al. 
2003). 







/ 


A 


R 


Q 


73.1 


30.0* 


11.6 


2.0 


0.54 


0.67 



The core-knot separation (i?) is in mas, the precession velocity, 
f2 is in deg/d, and all other parameters are in degree. * 30.0 deg 
corresponds to the initial phase of the smallest ellipse in Fig 5, 
other initial phases from from to large are, O.Odeg, 20.0deg and 
15.0 deg respectively. 

Table 4. The parameters derived from the fitting of the 
trace of knots A and B. 





V 


/ 


A 


R 


1.0 


60.0 


72.5 


30.7 


2.0 



The core-knot separation {R) is in arcsec, and all other parame- 
ters are in degree. 
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Distance (cm) 




Figure 1. The L-C of three ejecta corresponding to Eq. ((2]l. The 
three peaks corresponds to three knots at different distances to 
the core, Cl,Ca and C4. 



Figure 3. The observed (filled squares corresponding to Britzen 
et al. 2010) and fitted (opened circles) evolution of core separation 
(mas), flux and position angle (deg) of one component, CI of 
1803, in the time interval between 1984-1996. The corresponding 
from top to bottom panel are 6.3, 0.33, and 90.9 respectively. 




Figure 2. A schematic figure demonstrates knots produced by 
activities of the cental engine. Knots X1-X3 are produced by one 
pulse, and Y1-Y3 another. The upper panel corresponds to knots 
with small discrepancy in precession phases; and the knots of the 
bottom panel differ significantly in precession phases. 
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Figure 4. The observed and fitted ridge lines of BL Lac (Stirling 
et al. 2003). Points with error bars labeled A1-A4, B1-B4 and Cl- 
C4 are observed in December 1999, April 2000 and April 2001, 
respectively. The ellipses are predicted by the non-radial model. 
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Figure 5. The observation and interpretation of knots A and 
B of SS433. Top panel is the zeroth-ordcr HETGS image taken 
on 2001 May 12. Middle panel is the sketch of the orientation 
of the east jet of SS 433, which interprets the observations of 
knots A and B in context of the slowly moving clouds energized 
erratically by a more powerful, faster, unseen flow. The top and 
middle panels are adopted from Migliari et al. (2005). The bottom 
panel is the sketch explaining the observations of knots A and B 
by the non-ballistic model. 



Figure 6. A schematic explanation of the very large numbers 
of observed outward motions as opposed the inward ones. The 
curved jet precesses along the dash-dotted direction, and the ob- 
served knot can be seen as composed of several sub-knots, as 
shown in the left panel. The precession of the curved jet of the 
left panel is equivalent to a fix jet and a motion of the observer's 
LOS as shown in the middle panel. The right panel gives the 
observed knot, hollowed circle, and the core, filled circle. 



